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Abstract 

The full spectrum and eigenfunctions of the quantum version of a nonlinear oscillator de- 
fined on an TV-dimensional space with nonconstant curvature arc rigorously found. Since 
the underlying curved space generates a position-dependent kinetic energy three different 
quantization prescriptions are worked out by imposing that the maximal superintegrabil- 
ity of the system has to be preserved after quantization. The relationships among these 
three Schrodinger problems are described in detail through appropriate similarity trans- 
formations. These three approaches are used to illustrate different features of the quan- 
tization problem on TV-dimensional curved spaces or, alternatively of position-dependent 
mass quantum Hamiltonians. This quantum oscillator is, to the best of our knowledge, 
the first example of a maximally supcrintcgrable quantum system on an TV-dimensional 
space with nonconstant curvature. 
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1 Introduction 



This paper is devoted to the study of the quantum mechanical version of ^-dimensional 
(ND) classical Hamiltonian systems of the type 

«(q, P) = T(q, P) + W(q) = + W(q), 

where q, p G M. N are conjugate coordinates and momenta with canonical Poisson bracket 
{lijPj} = The physical interpretation of these systems is two- fold. On one hand, as 
the Hamiltonian describing the motion of a particle on the NT) curved space defined by the 
(conformally flat) metric ds 2 = 7W(q)dq 2 and under the action of the potential W(q). On 
the other, as position-dependent mass systems on the iVD Euclidean space. 

Evidently, the crucial point for the definition of the corresponding Schrodinger problem 
is the consistent (under certain given criteria) definition of the quantum kinetic energy term, 
T{<\, p) — > T(q, p), since an obvious ordering ambiguity appears when the position and 
momenta operators are considered. In this paper we shall deal with a specific example: the 
motion of a particle on the so-called NT) Darboux III space [HE] given by A^(q) = 1 + Aq 2 
(with A > 0), and whose potential will be the intrinsic oscillator on such space defined by 

ui \ ^ 2 q 2 

^ (q) = 2(TTAq^ 

This choice is motivated by the fact that this system is the only known example of a 
maximally superintegrable classical Hamiltonian on an NT) space with nonconstant curva- 
ture and it can be interpreted as a A-deformation of the flat isotropic oscillator, which is 
recovered in the A — > limit. In fact, it is well-known that the quantum superintegrability of 
the NT) flat isotropic oscillator is useful in order to obtain its exact solution by making use 
of the superabundance of quantum integrals of the motion Therefore, it seems natural to 
quantize the Darboux III oscillator in such a way that the maximal superintegrability of the 
system will be manifestly preserved under quantization; moreover, presumably in this way 
this new NT) nonlinear oscillator could be fully solved by mimicking the standard procedure 
for the Euclidean oscillator. 

In the sequel we will show that this is indeed the case. In fact, the superintegrability 
constraint will be useful in order to analyse in detail several possible quantization prescriptions 
amenable for iVD spaces with nonconstant curvature, namely: 

• The so-called [5] direct "Schrodinger" quantization 

which was the one used in [6] for this system, since it preserves the maximal superinte- 
grability in a straightforward way due to the immediate quantum transcription of the 
(2N — 1) classical integrals of the motion. This property leads to a maximal degeneracy 
of the spectrum, which is exactly the same as in the quantum NT) flat oscillator. 
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The "Laplace-Beltrami" (LB) quantization, which makes use of the usual LB operator 
on curved spaces: 

t2 2 2 N i 

H LB = —Alb + " q . 2 , where A LB = V —d^g^dj. (1.2) 
2 2(1 + Aq 2 ) f-± yfg 



However, we shall show that this LB Hamiltonian cannot be transformed into (1.1) 
through a similarity transformation unless we include an additional quantum potential 
which is proportional to the scalar curvature of the underlying space. This similarity 



transformation guarantees that the spectra of (1.2) and (1.1) coincide, and provides 
the explicit form of the full set of quantum integrals for %lb- We stress that such a 
kind of quantum "geometric" potential is well-known in the literature (both in scalar 
field theories in General Relativity as well as in the context of the quantization prob- 
lem on generic Riemannian manifolds 13 |9]), and it is tantamount to replace the 
LB operator by the so-called "conformal Laplacian" [10]. We have to mention here 
that the connection between LB operators and scalar curvatures associated with two 
different, possibly conformally flat, Riemannian manifolds has been firstly pointed out 
in a pioneering paper by Paneitz in 1983 [11]. However, the neat connection of the 
latter with quantum superintegrability properties is here stated for the first time and 
we think that this result opens the path for a novel algebraic approach to the subject. 

Finally, a "position-dependent mass" (PDM) quantization, which is essential in many 
condensed matter problems (see for instance [HI [131 021 EE E] ) • Here we perform 
it by considering the symmetric prescription proposed in [T5], namely: 

A x 1 . 1 . w 2 q 2 ft 2 1 w 2 q 2 . . 

?Wq, P ) = -p_^p+ 2(1 + Aq 2 ) = ~ T (TTAq 2 ) 2(1 + Aq 2 ) ' (0) 

In this case we again find that, in order to get a similarity transformation leading to 



(1.1), another additional quantum potential has to be added to (1.3). In this way the 
spectrum is preserved and the full set of quantum integrals of the motion is explicitly 
obtained. 



We stress that although the three previous quantum Hamiltonians have different explicit 
expressions and, to some extent, interpretations, all of them are related through similarity 
transformations provided the additional quantum potential terms are considered. This, in 
turn, means that they share a common energy spectrum, but they have different wave func- 
tions. Moreover, one of the main objectives of this comparative analysis is to point out some 
generic (i.e. potential independent) features of the quantization problem on curved spaces, 
such as the relevance of the dimension N of the underlying manifold (the N = 2 case will be 
distinguished) as well as the particular properties of the nonconstant curvature cases. 

The paper is organized as follows. In the next Section, the classical Darboux III oscillator 
is revisited in order to provide the necessary classical background on this system, including 
the explicit description of its maximal superintegrability in terms of a curved Fradkin tensor. 
Section 3 is devoted to review the geometry of the underlying space, emphasizing the role 
of its nonconstant curvature and introducing a classical radial effective potential that will 
be useful in the quantum context. In Section 4 the three superintegrable quantizations of 
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the Darboux III oscillator are obtained, and the similarity transformations among them 
are fully described. Section 5 provides the three associated radial Schrodinger equations. 
The spectrum and eigenfunctions of the system are rigorously obtained in Section 6, thus 
completing the preliminary results given in [jjj. Finally, Section 7 includes several remarks 
and open problems. 



2 The Darboux III oscillator 

The iVD classical Hamiltonian system given by 

with real parameters n > and oj > 0, was proven in [3] to be maximally superintegrable 
(MS), since this Hamiltonian is endowed with the maximum possible number of 2N — 1 
functionally independent constants of motion. Hereafter we shall consider the equivalent 
Hamiltonian T~L defined by % = kH/2 with real parameter A = 1/k > 0: 

H(q,p) = T( q ,p) + «( q) = 5Ir ^ + 5I1 ^. (2.2) 

The kinetic energy T(q, p) can be interpreted as the one generating the geodesic motion 
of a particle with unit mass on a conformally flat space with metric and (nonconstant) scalar 
curvature given by 

ds =(1 + Aq)dq, R(q) = —A (1 + Aq 2 )3 ^ ' 

In fact, such a curved space is the NT) spherically symmetric generalization of the Darboux 
surface of type III [TH1 EH], which was constructed in (TJ [2]. On the other hand, the central 
potential U was proven in [21 13] to be an "intrinsic" oscillator potential on that Darboux 
space. 



Moreover, in spite of the very naive appearance of the classical Hamiltonian (2.2), which 
is nothing but the iVD isotropic harmonic oscillator system divided by its oscillator potential 
(plus a rather relevant constant here scaled to 1), it is worth mentioning that this system can 
also be considered in three other different (but related) frameworks: 

• For N = 3, % arises as a particular case of the so call multifold (or i^-fold) Kepler 3D 
Hamiltonians constructed in [5j [20] as generalizations of the MIC-Kepler and Taub- 
NUT systems [2H E21 [231 (211 EH [261 123 [211 [29]. In the notation of [20J, U can be 
recovered by setting v = 1/2, a = 1 and b = A. Notice that from our approach, the 



proper Kepler-Coulomb potential on the Darboux space (2.3) would be [2] 



y/l + Aq 2 i i f~2 

«Kc = a rn > |q| = VQ i ® e 



which is related with the curved oscillator potential through U oc U K q ■ 
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Again for N = 3, the potential hi can be obtained from the temporal part of the 
family of the so called Bertrand metrics |30|, 131] on (3 + l)D Lorentzian spacetimes with 
nonconstant curvature. Such (3 + 1)D free systems possess stable circular orbits and 
all of their bounded trajectories are periodic, so these are the natural generalization of 
the classical Bertrand's theorem [32] to spaces of nonconstant curvature. We recall that 
the MS property for of all 3D Bertrand Hamiltonians, which come from the (3 + 1)D 
Bertrand metrics, was recently proven in [33J (see also 



As we have already mentioned, T~L can alternatively be interpreted as describing a PDM 



system in which the conformal factor of the metric (2.3) is identified with the mass 
function, which in this case is parabolic: .M(q) = 1 + Aq . 



We point out that we have chosen to deal with (2.2) instead of (2.1) because in this way 



all the expressions that we shall present throughout the paper will have a smooth and well 
defined limit A — > that leads to the well-known results concerning the (flat) NT) isotropic 
harmonic oscillator with frequency uj. In particular, the limit A — > of (2.2) and (2.3) yields 



1 2,1 2 2 

2 P + 2^ q ' 



d S ^ 



dq 2 , 



R = 0. 



(2.4) 



2.1 Maximal superintegrability from a curved Fradkin tensor 

The fact that % is a MS Hamiltonian can be explicitly demonstrated as follows [U 135] . 



Theorem 1. (i) The Hamiltonian li. (2.2) is endowed with the following constants of motion. 



(2N — 3) angular momentum integrals: 

C (m) = {q . v . _ qm) 2 



where m 



l<i<j'<m 



, N and C^* = C\ 



N-m<i<j<N 



QjPi) 



(2.5) 



(JV)- 



N integrals which form the ND curved Fradkin tensor: 



hj =PiPj ~ (2XH(q,p) - u} 2 )qiqj, 



(2.6) 



where i, j = 1, . . . ,N and such that H. = \ Y^i=\ hi- 

(ii) Each of the three sets {H, C (m) }, {H, C (m) } (m = 2,...,N) and {I u } (i = 1, . . . , N) is 
formed by N functionally independent functions in involution. 

(Hi) The set {T~L, C^ m \ Ci m \, la} for m = 2, . . . , N with a fixed index i is constituted by 2N — 1 
functionally independent functions. 



Notice that the first set of 2N — 3 integrals (2.5) is the same for any central potential on 
any spherically symmetric space [2] since it is provided by an underlying sl(2, R) coalgebra 



symmetry (also by an so(iV)-symmetry), while the second one (2.6) comes from the specific 
oscillator potential that we consider here. The latter, in fact, correspond to a curved analog 
of the Fradkin tensor of integrals of motion [1] for the isotropic harmonic oscillator. We also 



recall that the Hamiltonian (2.2) together with both sets of integrals of (2.5) and (2.6) can 
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alternatively be obtained |35j from the free Euclidean motion by means of a Stackel transform 
or coupling constant metamorphosis (see [361 S3 ESI EH] arid references therein). 



Thus, in general, the latter integrals (2.6) do not exist for a generic central potential so 



that, in principle, the MS property is not ensured at all. From this viewpoint the iVD nonlin- 



ear oscillator Hamiltonian % (2.2) can be regarded as the "closest neighbour" of nonconstant 



curvature to the isotropic harmonic oscillator system (2.4) (with A = 0) as both share the 
same MS property. In fact, the real parameter A behaves as a "deformation" parameter 
governing the nonlinear behaviour of 7-1, and this parameter is deeply related to the variable 
curvature of the underlying Darboux space. 

2.2 Expressions in terms of hyperspherical coordinates in phase space 

The above results can also be expressed in terms of hyperspherical coordinates r, 9j , and 
canonical momenta p r ,P0j, (j = 1, ■ ■ ■ ,N —1). The N hyperspherical coordinates are formed 
by a radial- type one r = |q| G M + and N — 1 angles 9j such that Ok 6 [0, 2ir) for k < N — 1 
and 0jv-i G [0, vr). These are defined by 

J-l N-l 

qj = rcos6j sin^, 1 < j < N, q N = r sin0 fe , (2.7) 
fc=i k=l 

where hereafter any product YiT such that I > m is assumed to be equal to 1. The metric 



(2.3) now adopts the form 

ds 2 = (l + Ar 2 )(dr 2 + rW), (2. 
where d$7 2 is the metric on the unit (N — 1)D sphere S N ~ l 

N-l j-l 

dn 2 = de 2 ]Jsm 2 9 k . 

j=i k=i 

The relations between p and p T ,pe j read (1 < j < N) [2|: 
3~ l „„,, a i- 1 rri- 1 



cos Vi pQ t _ 1 pq . 



TT ■ n n , COS % ^ Y\{=l+l °k 

Pj = 11 Sin^feCOS^Pr. + — ^- 2^ ni _! — 

k=l 

AT-1 AT-1 t-iN-1 ■ n 

sin 6> fe p r + - > , -, cos 6i p 6l , (2.9) 

fc=l r i=i n m =ism0 m 



where from now on any sum suc h that I > m is assumed to be zero. From (2.9) we 
obtain that 

p 2 =p 2 + r~ 2 L 2 , (2.10) 
where L 2 is the total angular momentum given by 

N-l j-i 

fri t=i sm ° k 
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By introducing (2.7) and (2.9) in the Hamiltonian (2.2) we find 

T(r, Pr )+U{r). 



v 2 + r~ 2 L 2 
H{r, Pr ) = P :J. ; k + 



u; 2 r 2 



2(1 + Ar 2 ) 2(1 + Ar 2 ) 



(2.12) 



The integrals of motion Cr m \ (2.5) adopt a compact form (the remaining and ly have 



more cumbersome expressions): 

N-l 



3-1 



c 'm= e pi n m = 2 



j= iV-m+l £;=7V-m+l 



iV; 



and C(at) = L 2 , which is just the second-order Casimir of the so(iV)-symmetry algebra of a 
central potential. 

Furthermore, the complete integrability determined by the set of N functions { r hL,Ci m )} 
(m = 2, . . . , N) leads to a separable set of N equations, since each of them depends on a 
unique pair of canonical variables. These are the N — 1 angular equations 



a 



(fc-i) 



C(fc) (9 N -k+l , P8 N _ k+1 ) = Pe N _ k+1 + . 2 a 

sm V N _ k+ i 

C {N) (0 1 , Pei )=pl + C ^^L 2 , 
sin ti\ 



k = 3,...,N -1, 



(2.13) 



together with the single radial equation corresponding to the ID Hamiltonian (2.12). 



3 The Darboux space and the classical effective potential 



The underlying manifold of the classical Hamiltonian (2.2) is the iVD Darboux space with 



metric (2.3), whose kinetic energy corresponds to the geodesic motion on the complete Rie- 
mannian manifold Ai N - 



DiV 



g) with 



9ij 



(1 + Aq 2 ) Si 



(3-1) 



and provided that A > 0. The scalar curvature R(r) = -R(|q|) (2.3) coming from this metric 
is always a negative increasing function such that lim,.-^ R = and it has a minimum at 
the origin 

12(0) = -2XN(N- 1), 

which is exactly the scalar curvature of the iVD hyperbolic space with negative constant sec- 
tional curvature equal to — 2A (see figure[l]). Recall that the four Darboux surfaces are the only 
2D spaces of nonconstant curvature whose geodesic motion is (quadratically) MS, therefore 
they are the "closest" ones to the classical Riemannian spaces of constant curvature [181 [T9] . 



As far as the nonlinear radial oscillator potential U(r) (2.12) is concerned, we find that it 
is a positive increasing function of r, such that 



U(r) 



u 2 r 2 



2(1 + Ar 2 )' 



W(0) = 0, 



lim U(r) 

r— >oo 



2A' 



(3.2) 
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This potential is shown in figure [2] for several values of A. Consequently, in contrast with 
the (Euclidean) isotropic harmonic oscillator, U(r) yields a nonlinear behavior governed by 
A, which means that the oscillator potential has the asymptotic maximum w 2 /(2A). 

Nevertheless, since the underlying manifold Ai N is not flat, the interplay between the 
oscillator potential U{r) and the kinetic energy term is rather subtle. For this reason, the 
complete classical system can be better understood by introducing a classical effective po- 
tential. This can be achieved by applying the ID canonical transformation defined by 

P( r, P r) = -» , Q(r) = lrVT^ + ^M, (3.3) 
V 1 + Ar z ^ 2v A 



(where the new canonical variables fulfill {Q,P} = 1), to the radial Hamiltonian (2.12). 
Notice that Q{r) has a unique (continuously differentiable) inverse r(Q), on the whole positive 
semiline, that is, both r,Q G [0, oo) and dQ(r) = Vl + Ar 2 dr. In this way, we obtain that 

H(Q,P) = \p 2 + U eS (Q), U eS (Q(r)) = 2(1 + C ^ >2 + 2{ ff Xr2 y (3-4) 

where the constant cn > is the value of the integral of motion corresponding to the square of 
the total angular momentum C(jv) = L 2 (2.13). Hence the classical system can be described 



as a particle on a ID flat space under the effective potential U c s(Q(r)), which is represented 
in figure [3| 

The analysis of U c q shows that this is always positive and it has a minimum located at 
^min such that 



r 



Ac/v + J \ 2 c 2 N + tu 2 c N 

min — ^2 ' 



^off(Q(^min)) = — Ac/v + \J^ 2 C 2 N +LxJ 2 C N . (3.5) 
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Figure 2: The nonlinear oscillator potential (J3T2J) with u = 1 for A = {0, 0.02, 0.04, 0.06, 0.1} 
starting from the upper dashed line corresponding to the isotropic harmonic oscillator with 
A = 0. The limit r — > oo gives {+oo, 25, 12.5, 8.33, 5}, respectively. 




0.02, c N 
3.49 with W e ff(r mill ) 



100 
8.2 



Figure 3: The classical effective nonlinear oscillator potential (3.4) for A 
and u> = 1. The minimum of the potential is located at r m i n 
and U e ff(cc) = 25. The dashed line corresponds to the effective potential of the harmonic 
oscillator with A = with minimum V(es{ r mm) = 10 at r m \ n = 3.16. 
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Therefore, r m j n and U e ff(Q(r m in)) are, in this order, greater and smaller than those corre- 
sponding to the isotropic harmonic oscillator: 

A = 0: r min = ^> ^cff(Q(?"min)) = Uy/CN- (3.6) 

Moreover U e g has two representative limits: 

uj 2 

HmZ4ff(Q(r)) = +oo, lim U eS (Q(r)) = — , (3.7) 

r— ^0 r— >oo 2A 



the latter being the same of (3.2). Thus, this effective potential is hydrogen-like and one 
should expect that its quantum counterpart should have both bounded and unbounded states. 
The rest of the paper is devoted to solve such a quantum problem in full detail. 

4 Superintegrable quantizations of the Darboux III oscillator 

Let us consider the quantum position and momenta operators, q, p, with Lie brackets and 
differential representation given by 

d 

W,Pj] = ihSij, Qi = Qi, Pi = • (4- 1 ) 

oqi 

Hereafter we will use the standard notation 

a a \ . „ 2 a 2 d 2 

V= 7— , A = V 2 = 7 — + ••• + 



dqi'"' dq N J ' d 2 qi d 2 q N ' 

4.1 The Schrodinger quantization 

The so-called "direct" quantization approach can be summarized in the following way [6] 
(this result is worth to be compared with Theorem 1, but taking into account that hereafter 
the order of the terms becomes crucial). 



Theorem 2. Let H be the quantum Hamiltonian given by 

H= , l — rs p 2 + ^%^r = , 1 9 , (-ft 2 A + ^ 2 q 2 )- (4.2) 
2(1 + Aq 2 ) V 2(1 + Aq 2 ) 2(1 + Aq V 1 



Then: 

(i) H commutes with the following observables: 

• The (2N — 3) quantum angular momentum operators, 

C (m) = ^ {qiPj ~ qjPif , C(m)= Yl (&Pj-VPi) 2 ' ( 4 - 3 ) 

l<i<j<m N-m<i<j<N 

where m = 2, . . . , N and <j( N ) = Crm. 

• The N 2 operators defining the ND quantum Fradkin tensor, given by 

hj = PiPj ~ 2A<7igj7i(q, p) + UJ 2 qiqj, (4.4) 
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where i,j = 1, . . . , N and such that H = \ ^2iL\ hi- 

(ii) Each of the three sets {H, C*( m )} ; {H, C\ m) } (m = 2, . . . , N) and {4} (% = 1, • • • , N) is 
formed by N algebraically independent commuting observables. 

(Hi) The set C( m \ Cr m \,In} for m = 2, . . . , N with a fixed index i is formed by 2N — 1 
algebraically independent observables. 

(iv) Ti. is formally self-adjoint on the I? Hilbert space defined by the scalar product 



f(q)$(q)(l + Aq 2 )dq. 



(4.5) 



M 1 



Proof. Some points of this statement can be straightforwardly proven through the coalgebra 
symmetry |40| HTl |4~2] of the quantum Hamiltonian (4.2). Let us consider the sl(2,M) Lie 
coalgebra in the basis { J±, J3} with commutation rules, Casimir invariant and (nondeformed) 
coproduct given by 



[J3, ^4 



2i^J 4 



C 



[J3,J- 
1 

2 

A(J,) = Ji ® 1 + 1 
An TV-particle realization of sl(2, R) reads 



-2ih J_ 
(J+J_ + J_J + ) 



FJ-, J4 



- 7 2 

/ = +,-,3. 



4i/iJs 



J4 



q 2 , 



^3 = ^(q-P + P-q)=q-p - ^i^- 



Therefore, % (4.2) has an sl(2,K) coalgebra symmetry since it can be written as 



H 



1 



2(1 + AJ_) 



2(1 + \J_ 



(4.6) 

(4.7) 
(4.8) 

(4.9) 
(4.10) 



Hence, by construction, % commutes with the (2A^ — 3) observables and C( m ) (m = 

2,...,N) (|4.3D which come from the "left" and "right" m-th coproducts [HI 02] of the 
invariant (4.7), respectively, up to an additive constant h?m(m — 4)/4. Furthermore, the 
coalgebra approach also ensures that these are algebraically independent and that each set 
{H,C^ m '} and {ti : C( m )} is formed by N commuting observables (to be more precise, they 
are polynomially independent as operators in a Jordan algebra). 

Next, by direct computations it can be proven that the iV 2 observables (4.4) commute 
with T~L, and that the N (diagonal) observables la (i = 1, . . . ,N) commute amongst them- 
selves as well; it is obvious that the latter la are algebraically independent. Finally, it is 
also clear that any single la is algebraically independent with respect to the set of 2N — 2 
observables C( m \ CV TO )} (as it is when A = 0) □. 

We stress that, as a byproduct of the above proof, any quantum Hamiltonian defined as 
a function of (4.9), 



H = H(J+, J-,J 3 ) = H(P 2 , q 2 , q • P - iHN/2), 



(4.11) 



is endowed with the same st(2, M) coalgebra symmetry. This shows that this is quasi-MS [21 
|4"H I42j . that is, it commutes, at least, with the (2iV — 3) observables CJ( m ) and C( m y In this 
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respect, we remark that what makes the quantum Darboux III oscillator (4.10) very special 



is the existence of a quantum Fradkin tensor formed by the "additional" symmetries Iij . This 
algebraic property implies that the system is MS and, as we shall see, that its discrete energy 
spectrum is maximally degenerate. 



4.2 The Laplace— Beltrami quantization 

When dealing with curved spaces with metric and classical kinetic term given by 

N . N 

2 



N - iV 

ds 2 = 9ij(<l)dqi<iqj, T(q,p) = - ^ 9 K ' (^)PiPj 



i,3=l 



the LB operator 



A LB = E \ ll i) ^"-lfl"'l 

i,j=l v y 



can be used in order to define the quantum kinetic energy as 

h 2 



7LB(q,p) 



-A 



LB, 



where g 13 is the inverse of the metric tensor and g is its determinant (see, for instance, [1J 
43J). If we apply such LB quantization to the Hamiltonian (2.2) with metric tensor (3.1) we 
get 

^ 2 a w2 q 2 

^LB = - y A LB + 2(1 + Aq2) 

h 2 h 2 X(N - 2) 

2(1 + Aq 2 ) 2(1 + Aq 2 ) 2(q ) + 



oj 2 q 2 



2(1 + Aq 2 )' 



Then, % (4.2) and 7^lb only coincide in the case N = 2 (as it should be for any sperically 
symmetric space |39j ) and for N > 2 they differ by a momentum-dependent potential, namely: 



n L B = n + Ui, Wi(q,p) 



■ hX(N-2) 
■ 1 2(l + Aq^ (q ' P) ' 



where we have introduced the quantum variables (4.1). Notice that 1A\ is linear in H, so this 



of an additional symmetry of the type (4.4). 



term does not have any classical analog. This situation reminds what happens in the context 
of the so-called quasi- exactly solvable quantum m odel s [44!. On the other hand, although the 
Hamiltonian ?^lb commutes with the operators (4.3) (the quantum correction IA\ preserves 
the sl(2,R) coalgebra symmetry (4.6)-(4.9)), in this case there is no hint about the existence 



Nevertheless, it is possible to find a "superintegrable" LB quantization (in the sense that 
it does preserve the MS property) by adding a second potential term to li. (besides U\) thus 
conveying iV 2 additional integrals of the type (4.4) together with the separability property 
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in terms of the N "diagonal" ones. In order to achieve this result we will relate T~L and %lb 
through a similarity transformation. If we apply Hlb to the product exp(/(q)) 1 I'(q) we get: 

h 2 \(N - 2W , „ T , h 2 X(N - 2)e' , „ „ , uiq 2 , , 

The two terms depending on (V/ • V^) and (q • V\P) can be removed by setting 

o _ /v 

/(q) = ^ r ln(l + Aq 2 ), (4.12) 
which, in turn, means that there is a similarity transformation connecting % and "Hlb : 



/i 2 A(iV - 2) 
8(1 + Aq 



^ = - ^1^3 (2* + 3Aq 2 (iV - 2)) . (4.13) 



Notice that the multiplication operator e-^ defines a transformation mapping 

L\M N ) = (R N , (1 + Aq 2 ) dq) into L 2 (M N ) = L 2 (r N , (1 + Aq 2 ) iV / 2 dq) , 



which is the natural L 2 space defined by the Riemannian metric. We remark that either 
for N = 2 or when A = 0, /(q) = W 2 (q) = 0, and notice also that the central potential 
IA2 is a pure quantum term as it depends on H 2 . The latter result suggests to consider a 
"transformed-LB" Hamiltonian defined by 

n T LB = n LB + u 2 = n+u 1 + u 2 , (4.i4) 

which satisfies 

n TLB =e f ne- f . (4.15) 

Hence, as a direct consequence, all the symmetries of % give rise to those corresponding to 
%tlb: 

X TLB = e f Xe~ f , X = {C<? n \C {m) J ij }, [Htlb, X tlb ] = 0. (4.16) 

Therefore, by taking into account Theorem 2 and the equations ( |4. 15 ) and (4.16) we find 
that "Htlb is, in fact, a quantum MS Hamiltonian. 

Theorem 3. Let %tlb be the quantum Hamiltonian given by 

Utlb = 2(1 + Aq 2 ) P + 2(l + A q 2 )-'2(l + A q ^ (c| - p) 
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Then: 



(i) "Wtlb commutes with the same observables (4-3), that is, CpLB = and CrLB,(m) 

C( m ) ; as well as with the N 2 Fradkin operators given by 

N — 2 N 



-^TLB, 



PiPj — (N — 2) 



2(1 + Aq 2 



;i + Aq 2 ) 2 



(N - 2)h 2 X 



2A<Mj% T LB(q, p) + ui Mj, 



(4.18) 



2(1 + Aq 2 ) ^ 

with i, j = 1, . . . , N and such that 7^tlb = h YliLi Itlb,u- 

(ii) Each of the three sets {T^tlb , } , {T^tlb, Cm} (m = 2, . . . ,N) and {Itlb,u} = 
1, . . . ,iV^ is formed by N algebraically independent commuting observables. 
(Hi) The set {7iTLB; C^ m \ C7 m ), Itlb,u} f or m = 2, . . . , N with a fixed index i is formed by 
2N — 1 algebraically independent observables. 

(iv) T^tlb is formally self-adjoint on the space L 2 (Ai N ) associated with the underlying Dar- 
boux space, defined by 



TLB 



*(q) $(q) (1 + Xq 2 ) N / 2 dq. 



(4.19) 



M' 



Therefore, according to the above statement, T^tlb can be seen as the appropriate LB- 



quantization of the classical Hamiltonian (2.2) as it manifestly preserves the MS property. 



Such a quantization requires to add a linear momentum-dependent potential U\ (coming from 
the quantum kinetic energy) plus an "additional" central potential XAi (coming from the MS 
property) to the Hamilt onian %. Clearly, the eigenfunctions of %tlb can be read off from 
those of % by means of (4.15). 



Remark. Interestingly, the quantum correction IA2 (4.13) to the oscillator potential aris- 



under lying metric (2.3): 



ing from the similarity transformation (4.12) is proportional to the scalar curvature of the 



Ui 



h 2 (N 



R, 



8(N-1) 

which vanishes for any 2D space and in the (most frequently studied) case of spaces of constant 
curvature gives simply an additional constant. Therefore, the transformed-LB quantization 
prescription is equivalent to imposing that 

h 2 



%TLB 



-A c +W(q) 



that is, to asserting that the appropriate quantum kinetic energy operator is essentially the 
conformal Laplacian (see, for instance, |1U| ) 

(N-2) 



A r =A 



LB 



4(N-1 



R, 



rather than the ordinary Laplacian (or LB operator). This is in full agreement with many 
prescriptions used in the analysis of scalar field theories in General Relativity or when dealing 
with quantization on arbitrary Riemannian manifolds However, we are not aware 

of any other instances where the convenience of this prescription has been motivated by 
superintegrability arguments. 
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4.3 A position-dependent mass quantization 

In the framework of PDM Hamiltonian systems there are also several ways to define the 
quantum kinetic energy term. A general approach depending on three parameters subjected 
to a constraint can be found in p3] (see also [E]). We shall consider here the proposal given 
in [15] and based on Galilean invariance arguments, which is the one extensively used in the 
condensed matter literature [T2], [13] . Such a PDM quantization, T(q, p) — > 7pDM(q>p), is 
defined as 2 

r PDM (q, p) = \ p ■ (rT \ ¥ yP = -y V • (j^sy V. 

Then, by adding the oscillator potential and ordering terms in the kinetic term, we obtain 
the following PDM quantization of the Hamiltonian (2.2): 

T^PDM = Tpdm (q, p) + U (q) 

h 2 h 2 X . „, u 2 q 2 



: A + m , ^ 2 (q-v) + 



2(1 + Aq 2 ) (l + Aq 2 ) 2VM y 2(1 + Aq 2 ) ' 



Hence the difference between H (4.2) and T^pdm relies again in a momentum-dependent 
potential, 

Wpdm = W + Vi, Vi(q, p) = i - A , 2 ^ 2 (q-p). 

Similarly to the LB quantization, the MS property can be explicitly restored through a 
similarity transformation and this process will require to add another central potential to the 
initial %pdm- 

Explicitly, if we apply T^pdm to the product exp(u(q))^(q) and define 

«(q) = ^ln(l + Aq 2 ), 

then we get the following similarity transformation between 7i and %pdm : 

H PDU e v V = e v (H- V 2 (q))*, ^pdm = e^e"" - V 2 , 

V2(q)= 2(lfA A q 2 )3 (^ + V(iV-3)). 

Hence, in contrast with the LB quantization, now both u(q) and V 2 (q) are nontrivial for any 
dimension N (including N = 2). In this way, we define the following "transformed-PDM" 
Hamiltonian, 

^TPDM=^PDM + V 2 = ?i + Vi+V 2 , Htpdm = e v He~ v , (4.20) 
whose symmetries are thus obtained from those of % as 

^tpdm = e^Xe"", X = {C^ m \ C( m ), I^}, [%tpdm, -^tpdm] = 0. (4-21) 
The MS property of the Hamiltonian Htpdm is summarized in the following statement. 



15 



Theorem 4. Let %tpdm be the quantum Hamiltonian defined by 

A 1 . 2 ^ 2 q 2 i^A , A M 7i 2 A(iV + Aq 2 (7V-3)) 

^TPDM - 2(1 - A ^ 2) P + 2(1 + A q 2) + (1 + M 2 )2 (q-P)+ 2(1 + Aq2)3 " 

(4.22) 

Then: 



(i) %tpdm commutes with the observables (4-3) as well as with (i,j = 1, . . . ,iV) 



-2A%g J ^ T PDM(q, p) + w 2 <Mj, 

which form a quantum Fradkin tensor and verifiy that %tpdm = \ Sj=i ^tpdm,h- 

(ii) Each of the three sets {%tpdm, C^}, {%tpdm, G(m)} (m = 2,...,N) and {/tpdm,m} 
(i = 1, . . . , N ) is formed by N algebraically independent commuting observables. 

(Hi) The set {T^tpdm, C^ m \ Cr m \ , Jtpdm,u} f or m = 2, . . . , N with a fixed index i is formed 
by 2N — 1 algebraically independent observables. 

(iv) T^tpdm is formally self-adjoint on the standard L? space with product 



(^|$)tpdm=/ #(q)$(q)dq. 

JM N 



Finally, we remark that by combining the similarity transformations (4.15) and (4.20) we 
obtain the relationship between %tlb and T^tpdm: 

^tpdm = e v -fn T LBe- {v - f) = (1 + Aq 2 ) Af / 4 7i T LB(l + \q 2 T m - 

5 Radial Schrodinger equations 

In this section we obtain the ID radial Schrodinger equation coming from each of the above 
three iVD quantum Hamiltonians by, firstly, i ntro ducing hyperspherical coordinates and, 
secondly, by making use of the observables Ct m \ (4.3 ) that encode the full spherical symmetry 
of the three systems. 



Let us introduce the map from the initial quantum operators (4.1 ) to the quantum hyper- 
spherical ones r, 0j, p r , pQ. (J = 1 , N— 1) with Lie brackets and differential representation 

given by 

[r,p r ] = iK, [frfOj] = 0, [Bj,Pr] = 0, [9j,P6 k ] = itrfjk, 

d - d 
f = r, p r = -ih—, 0j = d h p 6k = - 1 h—. (5.1) 

Here we point out that the "radial and phase operators" that we have just introduced are 
nothing but formal multiplicative operators on the angular variables, whose "canonical" trans- 
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formation rules with respect to the Cartesian ones are: 



r cos ( 



N-l 

QN = rY[ sin 4, 



k=l 



k=l 



i | ,, 4 - cosgj : ^ ]lfc = I + i sin 4 

Pi = 11 Sm 6k C ° S ^ Pr + — — n'- 1 .infl 

i JV-1 ttJV-1 • a 

i . 1 llk=l + l Smd k n » 

Pn = | | smV k p r + - > , -=jzi — r^— cos ^ P^i 



fe=i 

JV-1 

= n 



fc=i 

Hence we obtain that 



1=1 n m =i sin ' 



COS Ol p dl 



sm ( 



^nLlsin^fe 



P8i 



-2 -2 

q = r , 



1 



~JV-1^ L 2 ~2 -J^" 1 ) 
^Prr Pr + — = Pr -ifr 



Pr + 



q • P = rp r 



(5.2) 



where L 2 is the square of the total quantum angular momentum given by 

JV-1 fj-1 



1 



(sin 9j ) 



JV- 



— fyi&iOj)"- 1 -' fig 



Notice that the expressions (5.2) provide a ID (radial) representation of the sl(2,M) Lie 
algebra (4.6) by introducing them in (4.9). 



The N — 1 commuting observables CV m ) (4.3 ) turn out to be (m = 2, . . . , N) 

JV-1 



j-1 



a 



(m) 



e n 

j=AT-m+l \fc=JV-m+ 



1 



1 



r peAsm9j) N 1 ^ 



with C*(jv) = L 2 . Thus we obtain a set of N — 1 angular equations (& = 3, . . . , N — 1): 



C(k)(0N-k+l,P9 N _ k+1 



C(N){6liP6i] 



P0 JV _ fc+1 (sin0iv-fc + i) fe % JV _ fc+1 + 



(sin ( 



'JV-fc+1 



\fc-2 



C(fc-i) 



(smt/i 



iJV-2 



pepsin 6»i ) 7V 2 j5 9l + 



c, 



JV-2^ , ^(AT-1) _ f 2 

sm 9\ 



sm v N _ k+ i 
(5.3) 



which are worth to be compared with (2.13). Therefore the quantum radial Hamiltonian 
corresponding to (4.2) is obtained in the form 

1 



H(r,p r 



2(1 + Ar 2 ) \f 



1 * ~/V— 1 <- ^ 2 ^2 



(5.4) 



Afte r reordering terms and introducing the differential operators (5.1) in the Hamiltonian 
(5.4) we arrive at the following Schrodinger equation, %^ = E^f, 



1 



2(1 + Ar 2 ) 



-h 2 di 



2 -'> 2 h2{N 1} d r + ~ + ) tt(r, 0) = £¥(r, 0), 



(5.5) 
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where 6 = (9\, . . . , 9n-i). Next we factorize the wave function in the usual radial and 
angular components and consider the separability provided by the first integrals C( m ) (5.3) 
with eigenvalue equations given by 

*(r,0) = $(r)Y(0), (7 (m) * = c m *, m = 2,...,N. (5.6) 

Consequently, we obtain that Y(6) solves completely the angular part and such hyperspherical 
harmonics verify 

C {N) Y(0) = t 2 Y{0) = h 2 l(l + N -2)Y(9), 1 = 0,1,2... (5.7) 



where I is the angular quantum number. By taking into account the angular equations ( |5.3[ ), 
we find that the eigenvalues c m of the operators CV m ) are related to the N — 1 quantum 
numbers of the angular observables as 

Cfc -H- l k -i, k = 2, ...,N- 1, cjv-H-Z, 

that is, 

Hence the radial Schrodinger equation provided by % is 

1 / ^ ( d 2 (iV-l)d l(l + N-2)\ 9 9 \ _ . , . . 



2(1 + Ar 2 ) \ \dr 2 r dr r 2 , , 

In the same way, the ID radial Hamiltonian operators coming from the transformed LB 



(4.17) and PDM (4.22) quantizations are found to be 

~ h 2 / d 2 (N-l X(N-2)r\ d 1(1 + N -2) 

TLB ~ ~2(1 + Ar 2 ) VdT 2 + [~^~ + 1 + Ar 2 ) d7 r 2 

ft 2 / d 2 /JV-1 2Ar \ d 1(1 + N - 2) 



2(1 + Ar 2 ) \dr 2 \ r 1 + Ar 2 / dr r 2 

,2„2 



fr 2 A(iV + Ar 2 (iV-3)) 
2(1 + Ar 2 ) 2(1 + Ar 2 )3 ' W U; 

Recall that the three radial Hamiltonians %, T^tlb and %tpdMj are related through the 
similarity transformations as 

H TLB = (1 + Ar 2 )( 2 - 7V )/ 4 7i(l + Ar 2 )^" 2 )/ 4 , 
7i TPDM = (l + Ar 2 ) 1 / 2 ?i(l + Ar 2 r 1 / 2 , 

^tpdm = (1 + \r 2 ) N / 4 H TLB (l + \r 2 )~ N /\ 

Therefore the three corresponding radial Schrodinger equations share the same energy spec- 
trum and have different but equivalent radial wave functions: 

n<f>(r) = E$(r), 7i T LB$TLB(r) = E$ TLB (r), #TPDM$TPDM(r) = E<f> TPBU (r), 

$TLB(r) = (1 + Ar 2 )( 2 " Ar )/ 4 d>(r), $ TPDM (r) = (1 + Ar 2 ) 4 / 2 <I>(r), 

^TPDM(r) = (1 + Xr 2 ) N ^ ThB (r). (5.11) 
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6 Spectrum and eigenfunctions 



In this section we shall compute, in a rigorous manner, the (continuous and discrete) spectrum 
and ei genfun ctions of the quantum nonlinear oscillator by using the quantum Hamiltonian 



%tlb (4.22) characterized in Th eorem 3. We recall that the results corresponding to the 



Schrodinger quantization T~L (4.2) of Theorem 2 were advanced in [6] but without explicit 



proofs. Recall that both quantizations share the same spectrum but they have different 



radial wave functions which are related through the similarity transformation (5.11) 



6.1 Continuous spectrum 

Since M N is a complete manifold and the potential is continuous and bounded, it is standard 
that %tlb is essentially self-adjoint on the space Cq°(M n ) of smooth functions of compact 
support. It should be remarked that one cannot immediately determine the continuous 
spectrum of %tlb from asymptotics of the potential: in a complete Riemannian manifold, 
even the spectrum of the LB operator can be extremely difficult to analyze; e.g., it can be 
either purely continuous (as in Euclidean space), purely discrete [13] or consist of both a 
continuous part and eigenvalues, possibly embedded in the continuous spectrum [4"6] . 

In fact, to compute the continuous spectrum of %tlb it is convenient to take advantage 
of the spherical symmetry to decompose 

L 2 (M N ) = @L\R + ,dv)®y l , (6.1) 

zeN 

where dz/(r) = r N ~ 1 (l + Xr 2 ) N ' 2 dr and 3^z is the finite-dimensional space of (generalized) 
spherical harmonics, defined by 

y l ■= [Y G L 2 ^- 1 ) : A sN -iY = -1(1 + N- 2)Y) , 

where N stands for the set of nonnegative integers and Agjv-i denotes the Laplacian on the 
(N - 1)D sphere S^" 1 (or minus the angular momentum operator). This decomposition is 
tantamount to setting 

*TLB(q) = Yl Y ^°) $ TLB,zM , 

ZeN 

with = q/r£ S^ -1 , r = |q| and Y t £ y h 

As T^tlb is spherically symmetric, the decomposition (6.1) allows us to write Htlb as 
the direct sum of operators 

^TLB = 0£TLB,2®id y! , (6.2) 

ZeN 

with each Htlb,i standing for the Friedrichs extension of the differential operator on L 2 (M + , di^); 
namely 

. h 2 d JV _ 1 d h 2 X{N - 2)r d h 2 l(l + N - 2) 

2 -EM - ~ r JV-i(i + Ar 2 )d7 r &r~ (1 + Ar 2 ) 2 d7 + r 2 (l + Ar 2 ) 

1 + Ar 2 4(1 + Ar 2 ) 3 v v " 
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The continuous spectrum of Htlb is most easily dealt with using this decomposition. Indeed, 
from (6.2) it is apparent that 

spec(7i T LB) = (J spec(i?TLB,0 • 

To understand the spectrum of Htlb,i we proceed to compute and analyze its associated 
quantum effective potential it e fi,i- For this purpose we apply the same change of variable 
Q = Q( r ) (3.3) used in the classical case, together with a change of the radial wave function 
^ > TLB,K r ) ^ u (Q( r ))- We require that these transformations map the Schrodinger equation 
-f^TLB^TLB^ = -E^TLB.Z into 

-y|2 +U eSil (Q)\u(Q) = Eu(Q). (6.3) 
This is achieved by setting 

r (l-N)/2 

® TLB ' 1 W = (1 + Ar 2 )( 7V-l)/4 U W ( 6 - 4 ) 

in the radial Schrodinger equation, thus yielding 

,1 ,s 1 f H 2 (8(1 + \r 2 ) - 5) h 2 f 1 ri . N(N-A)\ 2 ,\ 

^ W = 2(T+A^j (4Wf + ^ + " " 2) + ) + " ^ J • 

(6.5) 

The behavior of £V e fr,i is rather similar to that of the classical effective potential (3.4) (see 
figure [i]), that is, U e sj is a positive function with a unique minimum, whose expression is 
rather cumbersome and which for the harmonic oscillator reduces to 

2 y/l(l + N-2) + (N-l)(N -^V4 

A U . ' min " ' 

to 

Weff,K^mm) = Hu>y/l(l + N - 2) + (N - l)(N - 3)/4. (6.6) 

Similarly to the classical system, the values of r m ; n and £4ff,z( r min) are r espe ctively greater 
and smaller than those corresponding to the quantum harmonic oscillator (6.6), but U e sj has 
the same asymptotic behaviour, namely, 

oj 2 

limZ4ff,j(r) = +00, lim ti eS i(r) = — . (6.7) 

We remark that there is a single exceptional particular case for I = and N = 2 for which 
U c fij reads 

^ (r) = 2(1 + Xr 2 ) { Ar 2 (l + Xr 2 ) 2 +W r J " 

Thus lim r _;.o ^4ff,z = — °° and linv-^oo ZVcff.z = w 2 / (2A), so W e fT,z has no minimum and can take 
both negative and positive values. 
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Figure 4: The quantum effective nonlinear oscillator potential (6.5) for N 
10 and h 



3, A = 0.02, 

I = 10 and H = ui = 1. The minimum of the potential is located, at T m ~ m — 3.59 with 
^4ff,z( r min) = 8.52 and £4^/(00) = 25. The dashed line corresponds to the quantum effective 
potential of the isotropic oscillator with A = with minimum W c ff,z(jmm) = 10.49 at r m ; n = 
3.24. 



Since we have just related the nonnegative, self-adjoint second-order differential operator 



on the half-line Htlb,i to (6.3), standard results in spectral theory J37J Theorem XIII. 7.66] 
ensure that the eigenvalues of i?TLB,z are contained in (0, E^) and its continuous spectrum 
is absolutely continuous and given by [-Eoo,oo), where we have set 

Eoo = lrni U cS j = — . 

Altogether, this guarantees that the continuous spectrum of Htlb is 

s P ec cont(^TLB) = [w 2 /( 2A ) 5 °°) , 

and that there are no embedded eigenvalues. 

6.2 Discrete spectrum and eigenfunctions 

Let us now compute the eigenvalues and eigenfunctions of T^tlb- To begin with, let us denote 
by ip n (q) the nth eigenfunction of the ID harmonic oscillator which satisfies 



2 v dg 2 y K ' V 2 
The explicit expression of ip n in terms of Hermite polynomials is 



■ifj n (q) = exp ( g 2 ) H n (J^q) , ( 



21 



up to a normalization constant. 

Due to the relationship between the Schrodinger and LB quantizations (4.15) we have 
that ^TLBlq) = (1 + Aq 2 )( 2_Ar )/ 4, I / (q) and the eigenvalue equation 

^TLB^TLB(q) = E^ TLB (q) 

can also be written as (see ( |4.2[ )) 

(-h 2 A + fi 2 q 2 )*(q) = 2£<f(q) , (6.9) 

where 

n = \Ju? - 2XE. (6.10) 

Since %tlb has no embedded eigenvalues (as shown in the previous subsection), one can 
safely assume that oj 2 — 2XE > 0. The condition ^tlb G L 2 (A4 N ) translates, according to 
( |4~5| ), as 

j |*(q)| 2 (l + Aq 2 )dq<oc; 

in particular, is square-integrable with respect to the Lebesgue measure. Therefore, by the 
standard theory of the harmonic oscillator, there must exist some n E N such that 

( N 

E = m[n-\ 

V 2 

Substituting the formula for f2, taking squares and isolating E, one readily finds that any 
eigenvalue of Htlb must be of the form 



E n = -Xfr n + 




(6.11) 

Conversely, one can prove that E n is an eigenvalue of %tlb for any n G N. This is easily seen 
by taking any partition (rii)f =l C N such that n\ + • • • + = n and noticing that, by (6.8) 
and do) ), 

N fo 

*TLB(q) = (l + Aq 2 )( 2 - Ar )/ 4 nexp{-/3V/2} J ff ni (/? (?i ), /3 = yj -, (6.12) 



i=l 

"2/ \aN 



is an L (Ai ) solution of the equation "Htlb^tlb = -En^TLB- 

Together with the result of the previous subsection, this proves the following 



Theorem 5. Let %tlb be the quantum Hamiltonian (4-22). Then: 

2 

(i) The continuous spectrum of %tlb is given by [7=^,00). Moreover, there are no embedded 
eigenvalues and its singular spectrum is empty. 

(ii) %tlb has an infinite number of eigenvalues, all of which are contained in (0, ^-). Their 

2 

only accumulation point is %r, that is, the bottom of the continuous spectrum. 

(Hi) All the eigenvalues o/'Htlb are of the form (6.11), and ^tlb is eigenfunction o/7^tlb 
with eigenvalue E n if and only if it is given by a linear combination of the functions (6.12) 
with n,i G N and n\ H V njv = n. 
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Figure 5: The discrete spectrum (6.11) for < n < 25, N = 3, h = u = 1 and 
A = {0,0.01,0.02,0.04} starting from the upper dot line corresponding to the isotropic 
harmonic oscillator with A = 0; in the same order, Eq = {1.5, 1.48, 1.46, 1.41} and 
Eoq = {oo,50,25,12.5}. 

Therefore the bound states of this system satisfy 

= lim E n = — , lim (E n+1 - E n ) = 0. 
Such a discrete spectrum is depicted in figure [5] for several values of A. 



7 Concluding remarks 

Summarizing, we have presented a novel exactly solvable quantum nonlinear oscillator in 
N dimensions, that can be understood as a simultaneous "analytic" A-deformation of both 
the usual isotropic oscillator potential and the underlying space on which the dynamics is 
defined. It turns out that if both sides of the Hamiltonian (the manifold and the potential) 
are appropriately modified, the curved quantum system preserves all the superintegrability 
properties of the Euclidean one, and its full solution can be explicitly obtained by making 
use of the curved analogues of Fradkin operators. It is worth stressing that such an explicit 
solution could be of interest from the physical viewpoint, since a parabolic effective- mass 
function has been proposed in |48t |4~9] in order to describe realistic quantum wells formed by 
semiconductor heterostructures. 

On the other hand, the quest for the preservation of superintegrability under quantization 
seems to be also a valuable guideline in order to clarify the properties of different possible 
quantization recipes on generic Riemannian manifolds. In this sense, we think that the 
connection here presented between the conformal Laplacian approach and the MS property 
is worth to be investigated through, for instance, the study of the quantization of other MS 
systems on spaces of nonconstant curvature that have been recently characterized in the 
context of the generalized Beltrami theorem |33j. 



23 



Finally, we recall that the real parameter A = 1/k was restricted in [3] to take a positive 
value. However, the MS of the classical Hamiltonian stated in Theorem 1 does hold for 
negative A as well. Nevertheless, the underlying space and the oscillator potential change 
dramatically when A < (see [35]), and the corresponding quantum problem is currently 
under investigation by making use of the techniques here presented. 
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